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THE STRAIN OF A NON-GRAVITATING SPHERE OF VARIABLE 

DENSITY* 

BY 

L. M. HOSKINS 

§ 1. Introductory. 

The mathematical problem of the strain of an elastic sphere is of special 
physical interest because of its relation to the question of the elastic yielding of 
the earth to disturbing forces. The cases of the problem which have been 
solved, however, involve assumptions which depart from what is known or 
thought probable in the case of the earth. The solutions of LamiS and Kelvin 
apply to a homogeneous body acted upon by bodily forces which are known 
functions of the coijrdinates. This case disagrees with the facts as regards the 
earth in two particulars : (1) The earth is certainly far from homogeneous in 
density and is very likely far from homogeneous in elastic properties. (2) The 
important bodily forces of self-gravitation cannot be expressed as functions of 
the coordinates in advance of the solution, since the changes in these forces 
caused by the strain are of the same order of importance as the external dis- 
turbing forces ; it is only by assuming incompressibility that it becomes possible 
to express the gravitational forces without the use of unknown functions. 

In a previous paper f was given a solution of the problem for a compressible 
sphere initially of uniform density, taking account of the changes in the gravi- 
tational forces caused by the strain. 

In the present paper will be found a solution for the case of a non-gravitating 
sphere of which the density is a given function of the distance from the center, 
the elastic moduli being assumed to be uniform in value throughout the body. 

By combining these two cases, and introducing the additional assumption that 
the elastic moduli are functions of the distance from the center, we should have 
perhaps the nearest possible approximation to a completely general statement of 
the problem presented by the actual earth. While there is no difficulty in 
forming the differential equations for this general problem in a manner analog- 
•ous to that employed in the two cases which have been solved, a solution of 
these general equations has not thus far been obtained. 

* Presented to the Society, San Francisco Section, February 26, 1910. 

fTransactions of the American Mathematical Society, vol. 11 (1910), p. 203. 
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The case treated in this paper appears, however, to possess some interest and 
importance, being a generalization in one important particular of the problem 
which formed the basis of the original estimates of the rigidity of the earth. 

§ 2. Statement of Prohlem. 

The problem here treated is the following : 

To determine the strain of an elastic sphere of which the density is a given 
function of the distance from the center, under the action of bodily forces hav- 
ing a potential, but free from stress on the bounding surface. The potential 
may be any function which, when expressed in polar coordinates, is developable 
as to 6 and ^ in a series of spherical surface harmonics. 

§ 3. Equations of Equilibrium of Volume-Element. 

Using polar co(5rdinates r, ^, ^, let u^, Ug, u^ be the elastic displacements, 
•nr^ , Wg, w^ the rotation-components, A the cubical dilatation, E^, E^, E^ the 
components of bodily force per unit mass, X + fyti , /x the bulk-modulus and 
rigidity-modulus, p the density. 

The three equations of equilibrium for an element of volume are the following : 

^ dA 2fi r d . . ^ dw.l „ 

(1) iX+2,)^-^^[^^i^,s^n0)-^\ + pE=O, 

dA 2a fd-BT d 1 

(2) (^+2;.)-^--.f^|^^-^^(...,sin^)J + pF, = 0, 

in which 

(^) 2.^« = ,:iiirg[^^-5r(""*^^'^^)]' 

IV d , du 1 

Let W denote the potential of the bodily forces, and assume that 
(8) W=^R,S,, 
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in wliich S. denotes a spherical surface harmonic of order i, and H. a function 
of r. Then 

dr •' 

(9) ^. = ^?^^N 



i^. = 2 



r dd 



r sin 6 dcf) 



A solution of equations (1), (2) and (3) may be formed by putting together 
solutions corresponding to the separate terms in the value of W. We there- 
fore confine the following discussion to the solution corresponding to a single 
term H. S, . 

We now assume 

(10) u. = uS., Ui, = v-^TTT, ti^ = —. — Ti-^rr, 
^ ^ '^ •' * 80 ' * sin ^ 5^ ' 

in which u and v are functions of r only. Substituting these values in equa- 
tions (4), • • • , (7), and making use of the differential equation 

(11) sin^^^'+^l(^-^W^)+^(^ + l)^^ = « 
which is satisfied by any spherical surface harmonic of order i, we find 

(13) 2.D-^ = 0,. 

« If d(rv)-\ 1 dS. 

(15) 2x^* = -4"~"^J"^- 

If the foregoing values of F^, Fg, F^, A, vr^, Wg and -ur^ be substituted in 
(1), (2) and (3), these are freed from and <p . Equations (2) and (3) in fact 
become identical, and there remain to be satisfied two ordinary differential equa- 
tions for determining m, z) as functions of r. Writing R for H^, these 
equations are the following: 



(16) 



„ , d fld(Pu) .,. ,,v\ 
('^ + 2m)^(;^-V-^(* + 1).-) 



i{i+\)t. (d{w) \.dR 
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/.rrx . o J'^d(r''u) .,. ^^v\ d fdirv) \ „ „ 

(17) ^X+2^.)[^-A^^ -^i^ + l)-y ^^^i^-^- -u) + pR = ^. 

It is now convenient to replace r hy x = r/a as independent variable, and to 
introduce y and z, defined by the equations 

.^o^ A 1 d(r^u) ., , _ ^v 

thus reducing the differential equations to the following : 

(20) (X + 2^)|+i.(z + l)M| + />g = 0, 

(21) (X + 2M)j/ + M^ + pi? = 0. 

§4. Surface Conditions. 

We now have to express the condition that the bounding surface is free from 
stress, that is, 

(22) rr = 0, rd = 0, r^ = (« = 1). 

Using (10) in the general formulas for these stresses, we find 

(23) ;;- = XA+2/~^ = (^X^ + 2/£'jS„ 

.... ->, /^w« Ug 1 5m \ fdv V u\dS. 

(24) ^^ = '^['af'--:^r'M)=^[d-r—r + 'ry8S^ 

/nr\ 'V / 1 ^^r ^'"'6 '"'d,\ ( dv V u\ 1 d 8. 

(2^) ''^ = ''\V^^e'd^+-dr-T) = ^\jr--r + -r)^e^'^ 
so that equations (22) reduce to the two 

(26) Xy + 2M§-0, 

to be satisfied when a; = 1 . 
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§ 5. Solution of Differential Equations. 

By eliminating successively s and y between equations (20) and (21) we obtain 
the two equations 

(28) (»,-|, + 24-i(i + l))((X+2M)y + ,i^)-^(.-Ji'|) = 0, 

(29) .(«-£ + 24-i(i + l)..)+««| = 0. 
The solutions of these equations are the following : 

( 2i + 1 ) [( X + 2/.)2/ 4- pR] = x* fx-'-' ^ L'E J) dx 

(30) 

— X-'-' i x'-T-( x'R-^) dx, 
J dx \ dx J 

(31) (2j + 1)ai.»= — x' j x-'-^ (xH -j-\ dx + x-*-^ I xMx-B^- Wx, 

in which, for the general solutions, each integral must be understood to involve 
an arbitrary constant. These constants are here omitted, since they are 
accounted for in the complementary solutions to be obtained presently. Modi- 
fying by integration by parts, we find 

(32) (X + 2^)y = - 2^1 [{i + l^fp- {§)dx + ^Jp^~^^dx], 

1 r . r d /^\ , 1 r d(Rx'+^) , -\ 

(33) ^. = 2i-^^[^a,J P^[^.)dx-^,j P-^dx'^''\- 

These values of y and z may now be substituted in (18) and (19), giving two 
simultaneous linear differential equations for determining u and v . The solu- 
tions of these equations may conveniently be expressed as follows : 

From (18) and (19) we readily obtain the equations 

dj i 4- 2 

(34) -^{iv-u)+-^{w-u) = a{iz-y), 

(^^^ ^[(^ + l)^ + ''*]-^[(^ + l)« + «] = «[(»■ + 1)^^ + 2/]' 

of which the solutions are 

(36) iv-u = ^,j{ iz - y)x'+'dx. 
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(37) {i + l)v + u = ax'-"^ i- 



'(^+i>±^cfe. 



These are easily solved for u and v . 

The forms of the complementary solutions are easily seen by considering the 
constants of integration in (32), (33), (36) and (37). The values of the con- 
stants must be determined so as to satisfy (18), (19), (20) and (21). The com- 
plementary solutions are thus found to be the following : 



(38) 



y = 


C.x* 


+ 


C-, 


-!«'-'■" 


-1 










z — 


B.x' 


+ 


-D-. 


_,x-- 


-1 










u 
r ~ 


A^' 


+ 


^-i- 


.,x-'- 


-' + A^_ 


_,x- 


-' + A_, 


-3«'-" 


-3 


r 


B.x* 


+ 


S-i- 


.,x-- 


-' + B, 


_,x'- 


-' + B_, 


.3X-- 


-3 



in which the twelve constants must satisfy the eight equations 

(\-h2/.)C, + (i + l)/.i?, = 0, 

A,-{i + 2)B, + D,= 0, 

(\ + 2,.)C_,_, -*/./)_,_, = 0, 

{i - 2)^_,_, + i{i + \)B_,_, + C_,_, = 0, 

A_,_, + {i-l)B_,_, + D_,_, = 0, 

A^_,-iB,_, = 0, 

A_,_, + ii + l)B_,_, = 0. 

The complete solutions wiU thus involve four arbitrary constants ; by as- 
signing values to these we may satisfy four equations of condition at the 
boundary. 

This solution covers the case of a spherical shell with bounding surfaces free 
from stress, or with arbitrary surface stress of a restricted type (corresponding 
to the restriction imposed on the potential of the bodily forces), or with arbitrary 
surface displacements of a similarly restricted type. 

For the case of the complete sphere, negative powers of x must be rejected, 
so that the only non-vanishing coefficients in (38) are A._2, A., B^_^, B., (7j, 
I)^ ; these being connected by four relations included among (39). There remain 
two arbitrary constants for satisfying the surface conditions (26) and (27). 

Trans. Am. Math. Soc. 33 
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§ 6. Case in which the Potential is a Harmonic Function. 

If TF is a harmonic function, each term H. S^ in its development must satisfy 
Laplace's equation, so that 

i?. = J.r' + Br-'-K 

Rejecting the negative power, we assume H. or S to have the form Ar'. 
Equations (32) and (33) then reduce to 

(40) (^■^^^)y = -(w^Jp--^P^-^ 

(^1) '^^=-(-2l+T)^J''-^.r-^^- 

The further discussion will be restricted to the ease in which i = 2 . 

§ 7. Sphere Strained hy Centrifugal Force. 

If the bodily forces are the centrifugal forces due to uniform rotation about a 
fixed axis, we have 

(42) ^=-|f-^A, 

in which c denotes the ratio of equatorial centrifugal force to gravity, and P^ is 
the zonal harmonic | cos^ ^ ~\- ^® ^^"^ t^^® 

(43) ,S,_-|P,_i-cos>«, 
,44) 5,„«_|^_^f^. 

In accordance with the notation of a former paper,* let e denote the ellipticity 
of a surface originally spherical and of radius r , and a the angular displax3ement 
of a radius vector for which ^ = 45° . With the above value of S^, we have 

(45) e = -, a=-. 

In writing the complete solutions of the differential equations for this case, 
we take A^ and A^ as arbitrary constants, expressing B^, B^, C^, D^ in terms 
of them by means of (39). We thus obtain the equations 

*Transactions American Mathematical Society, vol. 11 (1910), p. 206. 
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(47) ^ = ^^^3^A^-£3X%.^^., 

e-A +Ax^ + '^ r 5(3X + 5M) p 

(48) 

- ^^^^ Cpx'dx - 2(3X + Sfi) r pxdxl , 

« = i^o + — 6X — '^ 
(49) 

+ 30KX+2M) [Si f^'^''+ ^^^I P''^«'-(3^+8M)f P^^^] ^ 

Substituting these values in (26) and 27), we obtain the following equations for 
determining A^ and A^: 

2 A -A + 3o;^f^) [- 10 ( 9X + 10^) J' px^dx 
(50) 



+ 



72(X + yn) r px^dx~^{8\ + 8iJi.)C pxdxl = 0. 



A + 



8X+ Ifi 



-^^ + 30MX + 2M) [^(^^ + '^\f^"^^" 



«^ 3X 
(51) 

-24(X + yu,) r /w;«(;a;-2(3X+8/t) C pxdx'\=0. 

Solving, we find 

^0 = mr(xT2^)Jl9xTU^)[^^(^'^' + ^6^^^ + 28^^) £px^dx 
(52) 

-504(X + /t)-^ r/9a;«(?a; + 2(3X+8yit)(19X + 14yit) f pxdx\ 

(^^) ^= = 30M(x+2/IKl9'^^-l4M) -^^^^(^ + ^) (I ^'''^"' -J '^''^^) • 

Of especial interest are the surface values of e and a : 

(^^) '' ^ 2K1^X ^rn^T) [^(8^ + Tf^)£poo*dx - 21(X + m) jVa''^^^], 

(^^) «^ = 4M(19ni4M) [^^r^'^^ + ^^^^ + ^)f ^''^^- 
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If p is constant it is easy to show that the foregoing equations give the 
following well-known results : 

§ 8. Laplace^ s Law of Density. 

To illustrate the application of the above formulas to a case of variable density, 
assume Laplace's law, expressed by the formula 

,.„ _ sin KX p^K^ svcucx 

^ ' " "" «x 3 ( sin K — K cos k) a:x ' 

in which p^ is the mean density and p^ the value of p at the center. We then 
find 

C^ p K 

(59) ( pxdx = w-r^- "^ r(l — cos«x), 

^ '' J^ '^ 3(sm «— /ccosk)^ ^ 

..s 1 r ,, pic r6-«V 6-K^x^. 1 

(60 ) —, I px* ax = s-r-i ^ — 2—2 — cos kx s— i^ — sm kx , 

^ Jo o(sm «: — /iccos«:) |_ k x^ K^xr J 

1 r ,, pic V120-&0ic'x^ + bic'x' . 

(61) -. px^dx = ^s-r-- — ^ N 5-y sir 

^ Jo 3(sm « — « cos*) L K^x^ 

120 - 20K2a;^ + k^cc* 



sin «x 



cos «a; 



]■ 

If the ratio of mean density to surface density is assumed to be 2 , we find 
a; =2.461, and 

(62) r /9x*cZa; = .168/3„, ^ px^dx:=.10?>p^. 

'Jo "0 

These values substituted in (54) and (55) give 

e^ 4.45\ + 3.61/ii p^ga 



(63) 



(64) 



c 2(19X + 14/i) p, ' 
a^ 2.35X+ 1.51/ti /J^g'a 



c ~4(l9X+14At) /u. 
Comparing with the surface values of e and a for a sphere of uniform density 
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^^ c ~2(19X + 14/i)' iM ' 

^ ^ c ~4(19X + 14/t)' /t ' 

it appears that the assumption of variable density decreases e^ by about 11 per 
cent and a^ by about 22 per cent. 

We may also compare the values at the center, which with Laplace's law of 
density are found to be 

/fi7N 2a„ e, A^ 124.4X^ + 3 49.1 X/. + 191.8/.' p^ga 

^^ c ~c~ c~ 30(X+2At)(i9X+14Ai) ' /^ ' 

while the corresponding values for uniform density p^ are 



(68) 



1a,_e^_ 4X + 3m p^ga 120X^ + 330X/. + 180/." p^ga 
c c 19X + 14/1 IX 30(X + 2/.)(19X+14At) ■ /t ' 



The central values of e and a are thus slightly increased by the assumption of 
variable density. 

§ 9. Gravitating Sphere. 

The results of the above solution cannot be applied directly to a gravitating 
body of size and mass comparable with those of the earth, since the changes in 
the gravitational forces caused by the strain will be of the same order of impor- 
tance as the disturbing forces. There is no difficulty in expressing the gravita- 
tional forces in terms of the strain and the density and thus forming the differ- 
ential equations which must replace (16) and (17) above given, but the solution 
of these equations appears to involve serious analytical difficulties.* 

*In expressing the equations of interior equilibrium in terms of the strain, for the case of a 
gravitating body, it is necessary to assume that the changes of stress and strain from their initial 
values are related in the same way as if the Initial condition were one of zero stress and strain. 
This assumption is the basis of the solution given for the case of a compressible sphere in the 
paper cited above, the reasoning being in this respect substantially identical with that previously 
applied to the case of incompressibility. (See A. E. H. Love, A Treatise on the Mathematical 
Tlieory of Elasticity, 2d ed., p. 249. ) It has sometimes been assumed that this reasoning is valid 
only in the case of incompressibility. Thus Love states that if the sphere is homogeneous and 
incompressible, the initial stress may be taken to be of the nature of a hydrostatic pressure, its 
value being then easily determined, and the terms which correspond to this initial stress in the 
equations of equilibrium being found to cancel with those representing the initial body-force. 
It seems clear, however, first that even in an incompressible body the initial condition may 
involve tangential stress, and secondly that the assumption of no initial tangential stress is not 
essential to the argument. Whatever values the initial stress-components may have, the terms 
representing them in the equations of equilibrium will cancel the terms representing the initial 
body-forces. 
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If the body is assumed incompressible and of uniform density, the effect of 
gravitation can be taken account of synthetically by a method employed by 
Kelvin.* This method consists in combining the results for the two cases of a 
gravitating fluid and a non-gravitating elastic solid. It does not appear that 
similar reasoning can be applied if either of the assumptions of incompressibility 
and uniform density is dropped, since the distribution of the forces of gravita- 
tion cannot be assumed to be the same for the strained elastic solid as for the 
fluid. 



The method of applying this principle in the solutions cited above, both for the case of incom- 
pressibility and for that of compressibility, consists in comparing the forces acting upon a 
volume-element in the strained condition with those initially acting upon the same volume-ele- 
ment. It would perhaps be more reasonable to deal with an individual element of matter, using 
the increments of the stress-components acting on that element in expressing the stress-strain 
relations by the usual formulas. The increment of stress for an individual element would be in 
equilibrium with the increment of body-force acting upon that element, which would consist in 
part of the change in the gravitational force due to the displacement of the attracted element. 
This would introduce additional terms into the equations of equilibrium, but would add no 
difficulty to their solution. 

* Thomson and Tait's Natural Philosophy, § 840, Ed. of 1890. 



